Extinct and extant termites reveal the fidelity of behavior fossilization in amber.

Fossils encompassing multiple individuals provide rare direct evidence of behavioral interactions among extinct organisms. However, the fossilization process can alter the spatial relationship between individuals and hinder behavioral reconstruction. Here, we report a Baltic amber inclusion preserving a female-male pair of the extinct termite species Electrotermes affinis. The head-to-abdomen contact in the fossilized pair resembles the tandem courtship behavior of extant termites, although their parallel body alignment differs from the linear alignment typical of tandem runs. To solve this inconsistency, we simulated the first stage of amber formation, the immobilization of captured organisms, by exposing living termite tandems to sticky surfaces. We found that the posture of the fossilized pair matches trapped tandems and differs from untrapped tandems. Thus, the fossilized pair likely is a tandem running pair, representing the direct evidence of the mating behavior of extinct termites. Furthermore, by comparing the postures of partners on a sticky surface and in the amber inclusion, we estimated that the male likely performed the leader role in the fossilized tandem. Our results demonstrate that past behavioral interactions can be reconstructed despite the spatial distortion of body poses during fossilization. Our taphonomic approach demonstrates how certain behaviors can be inferred from fossil occurrences.

Defensive glands in Stylotermitidae (Blattodea, Isoptera).

The large abundance of termites is partially achieved by their defensive abilities. Stylotermitidae represented by a single extant genus, Stylotermes, is a member of a termite group Neoisoptera that encompasses 83% of termite species and 94% of termite genera and is characterized by the presence of the frontal gland. Within Neoisoptera, Stylotermitidae represents a species-poor sister lineage of all other groups. We studied the structure of the frontal, labral and labial glands in soldiers and workers of Stylotermes faveolus, and the composition of the frontal gland secretion in S. faveolus and Stylotermes halumicus. We show that the frontal gland is a small active secretory organ in soldiers and workers. It produces a cocktail of monoterpenes in soldiers, and some of these monoterpenes and unidentified proteins in workers. The labral and labial glands are developed similarly to other termite species and contribute to defensive activities (labral in both castes, labial in soldiers) or to the production of digestive enzymes (labial in workers). Our results support the importance of the frontal gland in the evolution of Neoisoptera. Toxic, irritating and detectable monoterpenes play defensive and pheromonal functions and are likely critical novelties contributing to the ecological success of these termites.

Evidence of cospeciation between termites and their gut bacteria on a geological time scale.

Termites host diverse communities of gut microbes, including many bacterial lineages only found in this habitat. The bacteria endemic to termite guts are transmitted via two routes: a vertical route from parent colonies to daughter colonies and a horizontal route between colonies sometimes belonging to different termite species. The relative importance of both transmission routes in shaping the gut microbiota of termites remains unknown. Using bacterial marker genes derived from the gut metagenomes of 197 termites and one Cryptocercus cockroach, we show that bacteria endemic to termite guts are mostly transferred vertically. We identified 18 lineages of gut bacteria showing cophylogenetic patterns with termites over tens of millions of years. Horizontal transfer rates estimated for 16 bacterial lineages were within the range of those estimated for 15 mitochondrial genes, suggesting that horizontal transfers are uncommon and vertical transfers are the dominant transmission route in these lineages. Some of these associations probably date back more than 150 million years and are an order of magnitude older than the cophylogenetic patterns between mammalian hosts and their gut bacteria. Our results suggest that termites have cospeciated with their gut bacteria since first appearing in the geological record.

Molecular Phylogeny Reveals the Past Transoceanic Voyages of Drywood Termites (Isoptera, Kalotermitidae).

Termites are major decomposers in terrestrial ecosystems and the second most diverse lineage of social insects. The Kalotermitidae form the second-largest termite family and are distributed across tropical and subtropical ecosystems, where they typically live in small colonies confined to single wood items inhabited by individuals with no foraging abilities. How the Kalotermitidae have acquired their global distribution patterns remains unresolved. Similarly, it is unclear whether foraging is ancestral to Kalotermitidae or was secondarily acquired in a few species. These questions can be addressed in a phylogenetic framework. We inferred time-calibrated phylogenetic trees of Kalotermitidae using mitochondrial genomes of ∼120 species, about 27% of kalotermitid diversity, including representatives of 21 of the 23 kalotermitid genera. Our mitochondrial genome phylogenetic trees were corroborated by phylogenies inferred from nuclear ultraconserved elements derived from a subset of 28 species. We found that extant kalotermitids shared a common ancestor 84 Ma (75-93 Ma 95% highest posterior density), indicating that a few disjunctions among early-diverging kalotermitid lineages may predate Gondwana breakup. However, most of the ∼40 disjunctions among biogeographic realms were dated at <50 Ma, indicating that transoceanic dispersals, and more recently human-mediated dispersals, have been the major drivers of the global distribution of Kalotermitidae. Our phylogeny also revealed that the capacity to forage is often found in early-diverging kalotermitid lineages, implying the ancestors of Kalotermitidae were able to forage among multiple wood pieces. Our phylogenetic estimates provide a platform for critical taxonomic revision and future comparative analyses of Kalotermitidae.

The functional evolution of termite gut microbiota.

BACKGROUND: Termites primarily feed on lignocellulose or soil in association with specific gut microbes. The functioning of the termite gut microbiota is partly understood in a handful of wood-feeding pest species but remains largely unknown in other taxa. We intend to fill this gap and provide a global understanding of the functional evolution of termite gut microbiota. RESULTS: We sequenced the gut metagenomes of 145 samples representative of the termite diversity. We show that the prokaryotic fraction of the gut microbiota of all termites possesses similar genes for carbohydrate and nitrogen metabolisms, in proportions varying with termite phylogenetic position and diet. The presence of a conserved set of gut prokaryotic genes implies that essential nutritional functions were present in the ancestor of modern termites. Furthermore, the abundance of these genes largely correlated with the host phylogeny. Finally, we found that the adaptation to a diet of soil by some termite lineages was accompanied by a change in the stoichiometry of genes involved in important nutritional functions rather than by the acquisition of new genes and pathways. CONCLUSIONS: Our results reveal that the composition and function of termite gut prokaryotic communities have been remarkably conserved since termites first appeared ~ 150 million years ago. Therefore, the "world's smallest bioreactor" has been operating as a multipartite symbiosis composed of termites, archaea, bacteria, and cellulolytic flagellates since its inception. Video Abstract.

The ultrastructure of the intramandibular gland in soldiers of the termite Machadotermes rigidus (Blattodea: Termitidae: Apicotermitinae).

Machadotermes is one of the basal Apicotermitinae genera, living in tropical West Africa. Old observations suggested the presence of a new gland, the intramandibular gland, in Machadotermes soldiers. Here, by combining micro-computed tomography, optical and electron microscopy, we showed that the gland exists in Machadotermes soldiers only as an active exocrine organ, consisting of numerous class III cells (bicellular units made of secretory and canal cells), within which the secretion is produced in rough endoplasmic reticulum, and modified and stored in Golgi apparatus. The final secretion is released out from the body through epicuticular canals running through the mandible cuticle to the exterior. We also studied three other Apicotermitinae, Indotermes, Duplidentitermes, and Jugositermes, in which this gland is absent. We speculate that the secretion of this gland may be used as a general protectant or antimicrobial agent. In addition, we observed that the frontal gland, a specific defensive organ in termites, is absent in Machadotermes soldiers while it is tiny in Indotermes soldiers and small in Duplidentitermes and Jugositermes soldiers. At last, we could also observe in all these species the labral, mandibular and labial glands, other exocrine glands present in all termite species studied so far.

Unmapped RNA Virus Diversity in Termites and their Symbionts.

Despite their ecological importance, nothing is known about the diversity and abundance of RNA viruses in termites (Termitoidae). We used a metatranscriptomics approach to determine the RNA virome structure of 50 diverse species of termite that differ in both phylogenetic position and colony composition. From these samples, we identified 67 novel RNA viruses, characterized their genomes, quantified their abundance and inferred their evolutionary history. These viruses were found within or similar to those from the Togaviridae, Iflaviridae, Polycipiviridae, Flaviviridae, Leviviridae, Narnaviridae, Mitoviridae, Lispivirdae, Phasmaviridae, Picobirnaviridae and Partitiviridae. However, all viruses identified were novel and divergent, exhibiting only 20% to 45% amino acid identity to previously identified viruses. Our analysis suggested that 17 of the viruses identified were termite-infecting, with the remainder likely associated with the termite microbiome or diet. Unclassified sobemo-like and bunya-like viruses dominated termite viromes, while most of the phylogenetic diversity was provided by the picobirna- and mitovirus-like viruses. Of note was the identification of a novel flavi-like virus most closely related to those found in marine vertebrates and invertebrates. Notably, the sampling procedure had the strongest association with virome composition, with greater RNA virome diversity in libraries prepared from whole termite bodies than those that only sampled heads.

Ocean currents promote rare species diversity in protists.

Oceans host communities of plankton composed of relatively few abundant species and many rare species. The number of rare protist species in these communities, as estimated in metagenomic studies, decays as a steep power law of their abundance. The ecological factors at the origin of this pattern remain elusive. We propose that chaotic advection by oceanic currents affects biodiversity patterns of rare species. To test this hypothesis, we introduce a spatially explicit coalescence model that reconstructs the species diversity of a sample of water. Our model predicts, in the presence of chaotic advection, a steeper power law decay of the species abundance distribution and a steeper increase of the number of observed species with sample size. A comparison of metagenomic studies of planktonic protist communities in oceans and in lakes quantitatively confirms our prediction. Our results support that oceanic currents positively affect the diversity of rare aquatic microbes.

Desaturase specificity is controlled by the physicochemical properties of a single amino acid residue in the substrate binding tunnel.

Membrane fatty acyl desaturases (mFAD) are ubiquitous enzymes in eukaryotes. They introduce double bonds into fatty acids (FAs), producing structurally diverse unsaturated FAs which serve as membrane lipid components or precursors of signaling molecules. The mechanisms controlling enzymatic specificity and selectivity of desaturation are, however, poorly understood. We found that the physicochemical properties, particularly side chain volume, of a single amino acid (aa) residue in insect mFADs (Lepidoptera: Bombyx mori and Manduca sexta) control the desaturation products. Molecular dynamics simulations of systems comprising wild-type or mutant mFADs with fatty acyl-CoA substrates revealed that the single aa substitution likely directs the outcome of the desaturation reaction by modulating the distance between substrate fatty acyl carbon atoms and active center metal ions. These findings, as well as our methodology combining mFAD mutational screening with molecular dynamics simulations, will facilitate prediction of desaturation products and facilitate engineering of mFADs for biotechnological applications.

Evolution of Termite Symbiosis Informed by Transcriptome-Based Phylogenies.

Termitidae comprises ∼80% of all termite species [1] that play dominant decomposer roles in tropical ecosystems [2, 3]. Two major events during termite evolution were the loss of cellulolytic gut protozoans in the ancestor of Termitidae and the subsequent gain in the termitid subfamily Macrotermitinae of fungal symbionts cultivated externally in "combs" constructed within the nest [4, 5]. How these symbiotic transitions occurred remains unresolved. Phylogenetic analyses of mitochondrial data previously suggested that Macrotermitinae is the earliest branching termitid lineage, followed soon after by Sphaerotermitinae [6], which cultivates bacterial symbionts on combs inside its nests [7]. This has led to the hypothesis that comb building was an important evolutionary step in the loss of gut protozoa in ancestral termitids [8]. We sequenced genomes and transcriptomes of 55 termite species and reconstructed phylogenetic trees from up to 4,065 orthologous genes of 68 species. We found strong support for a novel sister-group relationship between the bacterial comb-building Sphaerotermitinae and fungus comb-building Macrotermitinae. This key finding indicates that comb building is a derived trait within Termitidae and that the creation of a comb-like "external rumen" involving bacteria or fungi may not have driven the loss of protozoa from ancestral termitids, as previously hypothesized. Instead, associations with gut prokaryotic symbionts, combined with dietary shifts from wood to other plant-based substrates, may have played a more important role in this symbiotic transition. Our phylogenetic tree provides a platform for future studies of comparative termite evolution and the evolution of symbiosis in this taxon.

Expansion of the fatty acyl reductase gene family shaped pheromone communication in Hymenoptera.

Fatty acyl reductases (FARs) are involved in the biosynthesis of fatty alcohols that serve a range of biological roles. Insects typically harbor numerous FAR gene family members. While some FARs are involved in pheromone biosynthesis, the biological significance of the large number of FARs in insect genomes remains unclear.Using bumble bee (Bombini) FAR expression analysis and functional characterization, hymenopteran FAR gene tree reconstruction, and inspection of transposable elements (TEs) in the genomic environment of FARs, we uncovered a massive expansion of the FAR gene family in Hymenoptera, presumably facilitated by TEs. The expansion occurred in the common ancestor of bumble bees and stingless bees (Meliponini). We found that bumble bee FARs from the expanded FAR-A ortholog group contribute to the species-specific pheromone composition. Our results indicate that expansion and functional diversification of the FAR gene family played a key role in the evolution of pheromone communication in Hymenoptera.

Historical biogeography of the termite clade Rhinotermitinae (Blattodea: Isoptera).

Termites are the principal decomposers in tropical and subtropical ecosystems around the world. Time-calibrated molecular phylogenies show that some lineages of Neoisoptera diversified during the Oligocene and Miocene, and acquired their pantropical distribution through transoceanic dispersal events, probably by rafting in wood. In this paper, we intend to resolve the historical biogeography of one of the earliest branching lineages of Neoisoptera, the Rhinotermitinae. We used the mitochondrial genomes of 27 species of Rhinotermitinae to build two robust time-calibrated phylogenetic trees that we used to reconstruct the ancestral distribution of the group. Our analyses support the monophyly of Rhinotermitinae and all genera of Rhinotermitinae. Our molecular clock trees provided time estimations that diverged by up to 15.6 million years depending on whether or not 3rd codon positions were included. Rhinotermitinae arose 50.4-64.6 Ma (41.7-74.5 Ma 95% HPD). We detected four disjunctions among biogeographic realms, the earliest of which occurred 41.0-56.6 Ma (33.0-65.8 Ma 95% HPD), and the latest of which occurred 20.3-34.2 Ma (15.9-40.4 Ma 95% HPD). These results show that the Rhinotermitinae acquired their distribution through a combination of transoceanic dispersals and dispersals across land bridges.

Metabolomic and transcriptomic data on major metabolic/biosynthetic pathways in workers and soldiers of the termite Prorhinotermes simplex (Isoptera: Rhinotermitidae) and chemical synthesis of intermediates of defensive (E)-nitropentadec-1-ene biosynthesis.

Production of nitro compounds has only seldom been recorded in arthropods. The aliphatic nitroalkene (E)-nitropentadec-1-ene (NPD), identified in soldiers of the termite genus Prorhinotermes, was the first case documented in insects in early seventies. Yet, the biosynthetic origin of NPD has long remained unknown. We previously proposed that NPD arises through the condensation of amino acids glycine and/or l-serine with tetradecanoic acid along a biosynthetic pathway analogous to the formation of sphingolipids. Here, we provide a metabolomics and transcriptomic data of the Prorhinotermes simplex termite workers and soldiers. Data are related to NPD biosynthesis in P. simplex soldiers. Original metabolomics data were deposited in MetaboLights metabolomics database and are become publicly available after publishing the original article. Additionally, chemical synthesis of biosynthetic intermediates of NPD in nonlabeled and stable labeled forms are reported. Data extend our poor knowledge of arthropod metabolome and transcriptome and would be useful for comparative study in termites or other arthropods. The data were used for de-replication of NPD biosynthesis and published separately (Jirosová et al., 2017) [1].

Biotechnological potential of insect fatty acid-modifying enzymes.

There are more than one million described insect species. This species richness is reflected in the diversity of insect metabolic processes. In particular, biosynthesis of secondary metabolites, such as defensive compounds and chemical signals, encompasses an extraordinarily wide range of chemicals that are generally unparalleled among natural products from other organisms. Insect genomes, transcriptomes and proteomes thus offer a valuable resource for discovery of novel enzymes with potential for biotechnological applications. Here, we focus on fatty acid (FA) metabolism-related enzymes, notably the fatty acyl desaturases and fatty acyl reductases involved in the biosynthesis of FA-derived pheromones. Research on insect pheromone-biosynthetic enzymes, which exhibit diverse enzymatic properties, has the potential to broaden the understanding of enzyme specificity determinants and contribute to engineering of enzymes with desired properties for biotechnological production of FA derivatives. Additionally, the application of such pheromone-biosynthetic enzymes represents an environmentally friendly and economic alternative to the chemical synthesis of pheromones that are used in insect pest management strategies.

Co-option of the sphingolipid metabolism for the production of nitroalkene defensive chemicals in termite soldiers.

The aliphatic nitroalkene (E)-1-nitropentadec-1-ene (NPD), reported in early seventies in soldiers of the termite genus Prorhinotermes, was the first documented nitro compound produced by insects. Yet, its biosynthetic origin has long remained unknown. Here, we investigated in detail the biosynthesis of NPD in P. simplex soldiers. First, we track the dynamics in major metabolic pathways during soldier ontogeny, with emphasis on likely NPD precursors and intermediates. Second, we propose a hypothesis of NPD formation and verify its individual steps using in vivo incubations of putative precursors and intermediates. Third, we use a de novo assembled RNA-Seq profiles of workers and soldiers to identify putative enzymes underlying NPD formation. And fourth, we describe the caste- and age-specific expression dynamics of candidate initial genes of the proposed biosynthetic pathway. Our observations provide a strong support to the following biosynthetic scenario of NPD formation, representing an analogy of the sphingolipid pathway starting with the condensation of tetradecanoic acid with l-serine and leading to the formation of a C16 sphinganine. The C16 sphinganine is then oxidized at the terminal carbon to give rise to 2-amino-3-hydroxyhexadecanoic acid, further oxidized to 2-amino-3-oxohexadecanoic acid. Subsequent decarboxylation yields 1-aminopentadecan-2-one, which then proceeds through six-electron oxidation of the amino moiety to give rise to 1-nitropentadecan-2-one. Keto group reduction and hydroxyl moiety elimination lead to NPD. The proposed biosynthetic sequence has been constructed from age-related quantitative dynamics of individual intermediates and confirmed by the detection of labeled products downstream of the administered labeled intermediates. Comparative RNA-Seq analyses followed by qRT-PCR validation identified orthologs of serine palmitoyltransferase and 3-ketodihydrosphingosine reductase genes as highly expressed in the NPD production site, i.e. the frontal gland of soldiers. A dramatic onset of expression of the two genes in the first days of soldier's life coincides with the start of NPD biosynthesis, giving further support to the proposed biosynthetic hypothesis.

Regulation of Isoprenoid Pheromone Biosynthesis in Bumblebee Males.

Males of the closely related species Bombus terrestris and Bombus lucorum attract conspecific females by completely different marking pheromones. MP of B. terrestris and B. lucorum pheromones contain mainly isoprenoid (ISP) compounds and fatty acid derivatives, respectively. Here, we studied the regulation of ISP biosynthesis in both bumblebees. RNA-seq and qRT-PCR analyses indicated that acetoacetyl-CoA thiolase (AACT), 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR), and farnesyl diphosphate synthase (FPPS) transcripts are abundant in the B. terrestris labial gland. Maximal abundance of these transcripts correlated well with AACT enzymatic activity detected in the LG extracts. In contrast, transcript abundances of AACT, HMGR, and FPPS in B. lucorum were low, and AACT activity was not detected in LGs. These results suggest that transcriptional regulation plays a key role in the control of ISP biosynthetic gene expression and ISP pheromone biosynthesis in bumblebee males.

Molecular Mechanism of the Two-Component Suicidal Weapon of Neocapritermes taracua Old Workers.

In termites, as in many social insects, some individuals specialize in colony defense, developing diverse weaponry. As workers of the termite Neocapritermes taracua (Termitidae: Termitinae) age, their efficiency to perform general tasks decreases, while they accumulate defensive secretions and increase their readiness to fight. This defensive mechanism involves self-sacrifice through body rupture during which an enzyme, stored as blue crystals in dorsal pouches, converts precursors produced by the labial glands into highly toxic compounds. Here, we identify both components of this activated defense system and describe the molecular basis responsible for the toxicity of N. taracua worker autothysis. The blue crystals are formed almost exclusively by a specific protein named BP76. By matching N. taracua transcriptome databases with amino acid sequences, we identified BP76 to be a laccase. Following autothysis, the series of hydroquinone precursors produced by labial glands get mixed with BP76, resulting in the conversion of relatively harmless hydroquinones into toxic benzoquinone analogues. Neocapritermes taracua workers therefore rely on a two-component activated defense system, consisting of two separately stored secretions that can react only after suicidal body rupture, which produces a sticky and toxic cocktail harmful to opponents.

Evolution of moth sex pheromone composition by a single amino acid substitution in a fatty acid desaturase.

For sexual communication, moths primarily use blends of fatty acid derivatives containing one or more double bonds in various positions and configurations, called sex pheromones (SPs). To study the molecular basis of novel SP component (SPC) acquisition, we used the tobacco hornworm (Manduca sexta), which uses a blend of mono-, di-, and uncommon triunsaturated fatty acid (3UFA) derivatives as SP. We identified pheromone-biosynthetic fatty acid desaturases (FADs) MsexD3, MsexD5, and MsexD6 abundantly expressed in the M. sexta female pheromone gland. Their functional characterization and in vivo application of FAD substrates indicated that MsexD3 and MsexD5 biosynthesize 3UFAs via E/Z14 desaturation from diunsaturated fatty acids produced by previously characterized Z11-desaturase/conjugase MsexD2. Site-directed mutagenesis of sequentially highly similar MsexD3 and MsexD2 demonstrated that swapping of a single amino acid in the fatty acyl substrate binding tunnel introduces E/Z14-desaturase specificity to mutated MsexD2. Reconstruction of FAD gene phylogeny indicates that MsexD3 was recruited for biosynthesis of 3UFA SPCs in M. sexta lineage via gene duplication and neofunctionalization, whereas MsexD5 representing an alternative 3UFA-producing FAD has been acquired via activation of a presumably inactive ancestral MsexD5. Our results demonstrate that a change as small as a single amino acid substitution in a FAD enzyme might result in the acquisition of new SP compounds.

Δ12-Fatty acid desaturase from Candida parapsilosis is a multifunctional desaturase producing a range of polyunsaturated and hydroxylated fatty acids.

Numerous Δ12-, Δ15- and multifunctional membrane fatty acid desaturases (FADs) have been identified in fungi, revealing great variability in the enzymatic specificities of FADs involved in biosynthesis of polyunsaturated fatty acids (PUFAs). Here, we report gene isolation and characterization of novel Δ12/Δ15- and Δ15-FADs named CpFad2 and CpFad3, respectively, from the opportunistic pathogenic yeast Candida parapsilosis. Overexpression of CpFad3 in Saccharomyces cerevisiae strains supplemented with linoleic acid (Δ9,Δ12-18:2) and hexadecadienoic acid (Δ9,Δ12-16:2) leads to accumulation of Δ15-PUFAs, i.e., α-linolenic acid (Δ9,Δ12,Δ15-18:3) and hexadecatrienoic acid with an unusual terminal double bond (Δ9,Δ12,Δ15-16:3). CpFad2 produces a range of Δ12- and Δ15-PUFAs. The major products of CpFad2 are linoleic and hexadecadienoic acid (Δ9,Δ12-16:2), accompanied by α-linolenic acid and hexadecatrienoic acid (Δ9,Δ12,Δ15-16:3). Using GC/MS analysis of trimethylsilyl derivatives, we identified ricinoleic acid (12-hydroxy-9-octadecenoic acid) as an additional product of CpFad2. These results demonstrate that CpFAD2 is a multifunctional FAD and indicate that detailed analysis of fatty acid derivatives might uncover a range of enzymatic selectivities in other Δ12-FADs from budding yeasts (Ascomycota: Saccharomycotina).

The role of desaturases in the biosynthesis of marking pheromones in bumblebee males.

Bumblebee males (Hymenoptera) produce species-specific labial gland secretions called marking pheromones (MPs). MPs generally consist of terpenoids and fatty-acid-derived aliphatic compounds with various chain lengths predominantly containing one or no double bonds. The unsaturated fatty-acid-derived MP components were hypothesized to be produced by fatty acid desaturases (FADs) that exhibit diverse substrate specificities. To address this hypothesis, we isolated and functionally characterized FADs from three bumblebee species: Bombus lucorum, Bombus terrestris, and Bombus lapidarius. By employing RNA sequencing of the male labial glands and fat bodies of B. lucorum and B. terrestris, we identified five paralogous FAD-like sequences but only two FAD lineages were abundant and differentially expressed in the labial glands. We found that abundant FAD lineages were also expressed in the labial gland and fat body of Bombus lapidarius. Functional characterization of FADs in a yeast expression system confirmed that Δ4-FADs exhibited a unique Δ4-desaturase activity exclusively on 14-carbon fatty acyls and Δ9-FADs displayed Δ9-desaturase activity on 14- to 18-carbon fatty acyls. These results indicate that Δ9-FADs are involved in the biosynthesis of major unsaturated components of MPs in B. lucorum and B. lapidarius despite the diverse MP composition of these bumblebee species. The contribution of lipases, acyltransferases, esterases, and fatty acid reductases to production of the species-specific MP composition is also discussed in light of the transcriptomic data obtained in this study.